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ABSTRACT 

We report the detection of steady radio emission from the known X-ray source X9 in the glob¬ 
ular cluster 47 Tuc. With a double-peaked C IV emission line in its ultraviolet spectrum provid¬ 
ing a clear signature of accretion, this source had been previously classified as a cataclysmic 
variable. In deep ATCA imaging from 2010 and 2013, we identified a steady radio source at 
both 5.5 and 9.0GHz, with a radio spectral index (defined as cx ly^) of a = —0.4 zb 0.4. 
Our measured flux density of 42 zb 4 /iJy beam“^ at 5.5 GHz implies a radio luminosity {lyLy) 
of 5.8 X 10^^ergs significantly higher than any previous radio detection of an accreting 
white dwarf. Transitional millisecond pulsars, which have the highest radio-to-X-ray flux ra¬ 
tios among accreting neutron stars (still a factor of a few below accreting black holes at the 
same Lx), show distinctly different patterns of X-ray and radio variability than X9. When 
combined with archival X-ray measurements, our radio detection places 47 Tuc X9 very close 
to the radio/X-ray correlation for accreting black holes, and we explore the possibility that this 
source is instead a quiescent stellar-mass black hole X-ray binary. The nature of the donor star 
is uncertain; although the luminosity of the optical counterpart is consistent with a low-mass 
main sequence donor star, the mass transfer rate required to produce the high quiescent X-ray 
luminosity of 10^^ erg s“^ suggests the system may instead be ultracompact, with an orbital 
period of order 25 minutes. This is the fourth quiescent black hole candidate discovered to 
date in a Galactic globular cluster, and the only one with a confirmed accretion signature from 
its optical/ultraviolet spectrum. 

Key words: globular clusters: individual: 47 Tuc — X-rays: binaries — stars: black holes — 
radio continuum: general 


1 INTRODUCTION 

Globular clusters (GCs) are extremely efficient at producing low- 
mass X-ray binaries (LMXBs), with a specific frequency over two 
orders of magnitude higher than that seen in the field ( |Clark|1975l ). 
Although Galactic GCs comprise less than 10“^ of the total mass 
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of the Galaxy ([Harris, Harris & Alessi|2013[ [Licquia & Newmm 


2015 I, they contain roughly 10% of the LMXBs (Benacquista & 


Downing|2013|l. The high stellar densities of GCs imply additional 


formation pathways for LMXBs, via tidal capture and three body 
exchange interactions. However, while 18 dynamically-confirmed 
black holes (BHs) have been identified among the ~200 known 
LMXBs in our Galaxy (see [McClintock & Remillard|2006| for a 
review), none of these are in globular clusters. Either BHs are rare 
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in globular clusters, or selection effects have to date prevented their 
detection. 


1.1 Stellar-mass black holes in globular clusters 

The most massive stars in a newly-formed GC should collapse into 
BHs within the first few Myr, and it was originally thought that 
following mass segregation ( |Spitzer|1969^ , mutual gravitational in¬ 
teractions would lead to the expulsion of almost all of these BHs 


i|Kulkarni, Hut & McMillan|1993||Sigurdsson & Hernquist|1993| 

Portegies Zwart & McMillan 2000|l. However, recent detect! 

ons of 

in ex- 

variable X-ray sources in extragalactic GCs with luminosities 

cess of 10®® erg s“^ (e.g. Maccarone et al.|2007 Shih et al. 

2010 

Maccarone et al.|201 l[[Brassington et al.[2012[ Roberts et al. 

2012 


and references therein) have suggested that BHs can in fact survive 


in GCs. Closer to home, |Bamard et ^ ( |2014) have recently used 
X-ray spectral fitting to identify several BH candidates in GCs in 
M31. 

The implied survival of GC BHs has also received theoretical 
support (e.g. [Mackey et al.||2008| ), with recent simulations show¬ 
ing that GCs could in fact retain a significant population of BHs 
( [Sippel & Hurley|2013[|Morscher et al.|2013[ l. However, this pop¬ 
ulation of predicted BHs is sensitive to the magnitude of the natal 
kicks received by the BHs at formation. Assuming that they are not 
ejected at formation, the BHs will rapidly segregate to the cluster 
centre, and the exchange of energy between the BH suhcluster and 
the rest of the stars in the cluster allows a significant fraction of the 
BHs to be retained up to the present day. [Morscher et al.| ^2015^ 
performed a set of detailed Monte Carlo simulations, finding that a 
small fraction (~ 1 %) of the most massive black holes form a sub¬ 
cluster that undergoes a series of core oscillations, forming a steep 
cusp and ejecting the most massive single or binary BHs before re¬ 
expanding due to the interactions of three-body binaries with other 
objects. The rest of the (lower-mass BHs) remain permanently well 
mixed with the cluster stars. Depending on initial cluster parame¬ 
ters, up to 10^ BHs can remain after a Hubble time. However, once 
the BH population has been sufficiently depleted the BHs cannot 
provide sufficient energy to the rest of the cluster, which then en¬ 
ters a phase of “core collapse" (defined observationally via a steep 
radial density profile in the central regions), leading to the rapid 
ejection of the remaining BHs (|Breen & Heggie|2013| [Heggie &| 
|Giersz|2014> . This would imply that clusters that have not yet un¬ 
dergone this phase of observational core collapse should be the best 
candidates to host BHs. 

Although 18 LMXBs in GCs have been identified either as 
persistent sources or during transient outbursts (e.g. |Bahramian et| 
|al.|2014^ , sensitive X-ray observations have enabled the detection 
of dozens of candidate neutron star (NS) LMXBs in quiescence, 
identified by their soft thermal surface emission (e.g. [Rutle dg e et| 
[al.|2002[ [Heinke et al.|2003| [Guillot et al.|200^ , enabling the in¬ 
ference that there are of order 10-20 quiescent LMXBs for each 
LMXB seen in outburst, and thus > 200 quiescent NS LMXBs 
in GCs ( [Heinke et al.[[2005b^ . Identifying quiescent BH LMXBs 
in GCs is harder, since they do not produce a distinctive X-ray 
signature such as the thermal surface emission from quiescent NS 
LMXBs. However, accreting BHs do emit compact, partially self- 
absorbed radio jets in the hard and quiescent states (e.g. [Gallo,[ 
[Fender & H y nes[ [2005[(, with significantly higher radio luminosi¬ 
ties than NS systems with the same X-ray luminosity ( [Migliari'^ 
[Fender[[2006^ . Since the recent bandwidth upgrades to the Very 
Large Array (VLA) and the Australia Telescope Compact Array 
(ATCA) have now rendered this faint radio emission accessible 


over a significant fraction of the Galaxy, we can effectively search 
for BHs using deep radio continuum observations of GCs ( [Ma c-[ 
[carone[2004[[2005| (. The radio emission should be unresolved, with 
a flat or slightly inverted spectrum (i.e. a > 0, where we use 
the convention that flux density Su scales with frequency v as 
Sv oc i/“). The inferred radio luminosity should also fall on or 
within the scatter of the well-established correlation between X- 
ray and radio luminosities for hard state and quiescent accreting 
BHs jHannikainen et al.[1998[|^rbel et al.|2003|[Gallo, Fender &[ 

[ Pooley| 2003[(. The form of this correlation at the lowest luminosi¬ 
ties remains relatively poorly constrained ( [Gallo et al.|2006[ l, with 
recent evidence suggesting more scatter about the correlation than 
had previously been assumed ( [Gallo et al.|2014^ . However, in the 
absence of a dynamical mass, the presence of faint, flat-spectrum 
radio emission is still the best way to discriminate between an ac¬ 
creting BH and an analogous NS or white dwarf system. As origi¬ 
nally suggested by [Maccarone & Knigge[j2007[ l, this technique has 
recently been used to great effect, discovering two BH candidates in 


the core of M22 i 

Strader et al.|2012a|l and a third candidate in M62 

([Chomiuk et al.[ 

2013^. While follow-up observations to confirm 


the nature of these three flat-spectrum radio sources are ongoing, 
we have continued to obtain deep radio images of GCs to search 
for further BH candidates. 


1.2 47 Tucanae 


The dense, massive (6.5 x 10® Mq; Kimmig et al.|2015 i GC 47 
Tuc (NGC 104) lies at a distance of d = 4^57 ± 0.06 kpc (tak¬ 
ing the average of multiple previous distance measurements and 
the standard error on the mean; |Woodley et al.|2012^ . Owing to its 
proximity and low extinction {E{B — V) = 0.04 ± 0.02;|Salaris| 
[et al.|2007^ , the cluster has been very well studied, and is rich in 
exotic binary systems. 47 Tuc contains 300 known X-ray sources, 
including five quiescent NS LMXBs ( Heinke et al.|2005at and 23 
radio millisecond pulsars jCamilo et al. 12000 Freire et al.|2003[ l, 
which are likely LMXB descendants. Since conditions in the clus¬ 
ter are clearly favourable for forming LMXBs, this makes it an ideal 
target for deep radio observations to search for quiescent accreting 
BHs. Furthermore, [Giersz & Heggie[j201 1) recently performed de¬ 
tailed simulations of 47 Tuc, finding that it has not yet undergone 
core collapse (which will take another ~ 20 Gyr) and so should still 
retain a small population of BHs. This is in stark contrast to clusters 
such as M4 and NGC 6397 that have already undergone their “sec¬ 
ond core collapse" phase and would no longer be expected to host 
a significant BH population (see fig. 1 of [Heggie & Giersz|2014) . 

The first sensitive radio continuum observations of 47 Tuc 
were performed by [McConnell & Ables| j2000^ , who used 60 h 
of ATCA data at 1.4 and 1.7 GHz to reach rms noise levels of 
42 and 46/iJy beam“®, respectively. They detected eleven radio 
point sources within 5' of the cluster centre, including five sources 
later confirmed as millisecond pulsars (MSPs). By extending the 
on-sou rce integration time to 170 h, [McConnell, Deacon & Abies] 
(2001 I reached an rms noise level of 18 /rjy beam' , detecting an 


additional seven MSPs and three more unidentified sources, which 
they suggested to be either extragalactic background sources or 


as-yet unidentified MSPs. Within the 24 core radius (Howell, 


[Guhathakurta & Gilliland|||2000[l, the only detected sources were 
MSPs F, G/I, L and O. 

Following the recent bandwidth upgrade to the ATCA jWil-[ 
[son et aL][2011^ , [Lu & Kong[ ( [2011[ l made deep 5.5 and 9.0GHz 
observations of 47 Tuc in 2010 January, with a total on-source in¬ 
tegration time of 18.3 hr, aiming to detect the radiative signature of 
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a possible intermediate-mass black hole (IMBH) at the cluster cen¬ 
tre. After stacking both frequency bands to reach a noise level of 
13.3 /ijy heam“^, they detected no source brighter than So- within 
the central 10” x 10” of the cluster, giving a 3cr upper limit on the 
mass of any IMBH of 52O-49OOM0. A similar IMBH limit was 
found hy |McLaughlin et al.| ( |2006t , who used the cluster velocity 
dispersion profile to place a la upper limit of ISOOMq on the mass 
of any IMBH, so the evidence for such an object is not compelling. 

In this paper, we combine the archival observations of|Lu &| 
|Kong| j201 1^ with new ATCA data to make the deepest radio im¬ 
age of 47 Tuc to date. We report the discovery of a flat-spectrum 
radio source within the cluster core that is an astrometric match to 
the known bright X-ray source X9 (Hertz & Grindlay|1983[ |Au-| 


riere, Koch-Miramond & Ortolani|1989^ . In light of its inferred ra¬ 
dio luminosity, we suggest that the source is a quiescent BH, rather 
than a cataclysmic variable (CV), as had been previously inferred 
(|Paresce, de Marchi & Ferraro|1992||Grindlay et aL|2001||Knigge| 

let al.|2008ir 


In Sections and we describe our radio observations and 
results. We discuss previous observations of X9 in Section|^ com¬ 
paring the source properties to those of various classes of accreting 
compact object. We discuss the possible nature of the system in 
Section]^ review the sample of GCs hosting BH candidates in Sec¬ 
tion]^ and present our conclusions in Section|7] 


2 OBSERVATIONS AND DATA REDUCTION 

2.1 ATCA data 

We observed the globular cluster 47 Tuc with the ATCA on 2013 
November 12th, from 08:06-17:52 UT (MJD 56608.54 ± 0.20), 
achieving an on-source integration time of 8.7 hr. Using the Com¬ 
pact Array Broadband Backend (CABB), we observed simultane¬ 
ously in two bands, with central frequencies of 5.5 and 9.0 GHz, 
each with a bandwidth of 2048 MHz. The array was in its extended 
6A configuration, with a maximum baseline of 5.939 km. We used 
B1934-638 as both a bandpass calibrator and to set the flux density 
scale, and B2353-686 as the secondary calibrator to set the ampli¬ 
tude and phase gains. 

We reduced the data for each frequency band separately, 
performing external gain calibration in Miriad ( |Sault, Teuben| 
|& Wright||1995| l using standard procedures. We then frequency- 
averaged the calibrated data and imported them into the Common 
Astronomy Software Application (CASA; |McMullin et al.|2007^ 
for imaging. Imaging was carried out using Briggs weighting with 
a robust parameter of 1, which provided a good compromise be¬ 
tween sensitivity and resolution, as well as suppressing the side- 
lobes of the dirty beam. Our final image sensitivities of 5.7 and 
6.4/rjybeam“^ at 5.5 and 9.0GHz, respectively, were close to the 
theoretical thermal noise levels. The final source fluxes and posi¬ 
tions were measured in the image plane using the JMFIT algorithm 
within the Astronomical Image Processing System (AIPS; [Greisen| 
|2003l l. 

2.2 Archival ATCA data 

To improve the significance of our detections, we combined our 
data with the archival ATCA data taken in the same array configu¬ 
ration by |Lu & KongIpOl 1[ | on 2010 January 24-25 (summarised 
briefly in Section |1.2^ . (jur re-analysis of these data allowed us 
to reduce the noise levels in the 5.5 and 9.0 GHz images to the 
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Figure 1. 5.5-GHz ATCA image of the core of 47 Tuc. The red cross de¬ 
notes the cluster centre, the red circle shows the 24" core radius, and the 
thick orange circle highlights our detection of X9. 


theoretically-expected limits of 6.7 and 9.8/rJybeam~^, respec¬ 
tively. Despite the longer on-source integration times, these obser¬ 
vations were less sensitive than those from 2013 (Section [2. H be¬ 
cause the older data were taken prior to the installation of the new, 
more sensitive 4-cm receivers at the ATCA. 


3 RESUETS 

After stacking both the new and archival ATCA data in the 
uu-plane, our resulting images reached rms noise levels of 
4.4/iJybeam“^ at 5.5GHz and 5.7/ijybeam“^ at 9.0GHz. The 
final images at 5.5 and 9.0GHz had resolutions of 2'.'3 x l'.'7 in 
PA. 18°, and l"5 x I'.'l in PA. 28°, respectively. 


3.1 A radio counterpart to X9AV42 


The brightest radio source within the 24” core radius of the cluster 
at both frequencies (Fig.[^ is coincident with the second brightest 
hard X-ray source in the cluster (see Fig.®, denoted as source X9 
in the ROSAT catalogue of |Hasinger, Johnsto n & Verbunt (|1994| l, 
and W42 in the Chandra catalogues of |Griridl ay et al.| ( |2001| l and 
IHeinke et al.| ( |2005a[ >. We summarise in Table [T] the measured ra¬ 
dio brightness of the source in each observation. Using the stacked 
image, a point source fit to the radio source position in the image 
plane using the AIPS task JMFIT gave a J2000 position of 


R.A. = 00‘°24“04’;264 ± 0!016 
Dec. = -72°04“58”09 ± O'.'lO. 


Our fitted peak flux densities for X9 in the stacked data were 
42 ± 4 /ijy beam~^ at 5.5 GHz and 35 ± 6 /ijy beam“^ at 9.0 GHz. 
We therefore derive a radio spectral index of a = —0.4 ± 0.4 
and a 5.5-GHz radio luminosity of 5.8 x 10^^ ergs”'^. Here, and 
throughout this paper, we define radio luminosity by assuming a flat 
radio spectrum up to the observing frequency (i.e. 
where d is the source distance, v is the observing frequency, and Sv 
is the observed flux density). No circular polarisation was detected 
from the source in either of the individual epochs, with the best 5a 
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Figure 2. HST F814W image, showing the Icr error ellipses of 47 Tuc 
X9AV42/V1 in the radio (orange), optical (red) and X-ray (blue) bands. 
North is up and east is to the left. Uncertainty on the radio position arises 
purely from the statistical error in fitting the source position. Optical and X- 
ray error ellipses are dominated by the uncertainty in the absolute astrome¬ 
try (matching the HST reference frame to 2MASS in the optical band, and 
using the measured positions of X-ray counterparts to known radio MSPs 
to estimate the uncertainty on the X-ray reference frame). The en'or ellipses 
in all three wavebands overlap. 


Table 1. Measured flux densities of 47 Tuc X9. The source was fit as a point 
source in the image plane, using the AIPS task JMFIT. 


Data 

set 

(J2000) 

5.5-GHz 
flux density 
(/4Jy bm“^) 

9.0-GHz 
flux density 
(pjy bm“ 1) 

Spectral 

index 

2010 January 24-25 

34 ± 7 

36 ± 10 

0.2 ± 0.7 

2013 November 12 

50 ±6 

36 ±6 

-0.6 ±0.5 

Stacked data 

42 ±4 

35 ±6 

-0.4 ± 0.4 


limit in the stacked data set being < 22 /rjy beam“^ (< 52%) at 
5.5 GHz. We also placed a 5 (t limit of < 30 /rjy beam“^ (< 71%) 
on the 5.5-GHz linear polarisation. 

To quantify the variability of the radio counterpart to X9, we 
examined the individual images from the 2010 and 2013 data sets. 
Point source fits to the position of X9 showed the 9-GHz flux to 
be constant within uncertainties, while the 5.5-GHz flux increased 
from 34 ± 7 to 50 ± 6 /rjy beam“^, implying marginal (1.6a) ev¬ 
idence for variability and a slight steepening of the spectrum be¬ 
tween the two epochs (see Table[TJ. We also split the more sensitive 
2013 data set into two halves, finding no evidence for significant 
variability on timescales of a few hours. 


3.2 Other radio sources 


While X9 was the only source in the cluster core detected at > 5a, 
there were several other radio sources in the field. We used the 
AEGEAN source finder ( [Hancock et al.|2012a|bt to detect and char¬ 
acterise all sources within the half-power response of the primary 
beam at both frequencies (radii of 4!5 and 2'.1 at 5.5 and 9.0 GHz, 
respectively). AEGEAN fits elliptical Gaussian components (not 
constrained to the synthesised beam size, hence not necessarily 
point sources) to all islands of pixels with a peak exceeding a user- 
set threshold relative to the local background rms noise level, re¬ 
turning the amplitude of an elliptical Gaussian fit to each compo¬ 
nent. With over 10^ pixels in our images, we set a detection thresh¬ 
old of 5.5 (t. Including X9, AEGEAN detected fourteen radio sources 
exceeding this threshold within the central 4^5 FWHM of the 5.5- 
GHz image (Fig.[^. At 9.0GHz, only three sources (all with 5.5- 
GHz counterparts) were detected within the 2(7 diameter FWHM. 
Details of all detected sources are given in Table[^ 

We also examined the positions of the known MSPs in our 5.5- 
GHz image (identified from the ATNF pulsar catalogu^ |Manch- 


jester et al.|2005^ , and although none were detected at the 5.5a level 
(and hence picked up by AEGEAN), seven showed radio emission 
at > 3cr (Table [^. We carried out a point source fit to each, us¬ 
ing the JMFIT task within AIPS. For all other MSPs in the cluster, 
we place a 3(t upper limit of 15 /iJy beam“^ at 5.5 GHz. No MSPs 
were detected in the 9.0-GHz image. 

Of the 14 sources detected at 5.5 GHz by AEGEAN, four had 
an X-ray counterpart in the Chandra catalogue of|Grindlay et al.| 
( |200T) . These were W42 (X9), W132, W186 and W231, but only 
X9 had previously been identified with any particular source class. 
Given the flat radial distribution of X-ray sources beyond 100" 
(where the other three radio sources with X-ray counterparts are 
located), [Heinke et al.j ( |2005a[ l inferred that such distant sources 
are likely not associated with the cluster, but are more likely back¬ 
ground AGN. Indeed, optical imaging with the Hubble Space Tele¬ 
scope (HST) shows that W186 and W231 are likely associated with 
extragalactic sources since they since they lie in a region of ex¬ 
tended optical emission (van den Berg et al. in prep). Although 
W132 is outside the field of view of that data set, it is a very hard 
X-ray source ( [Heinke et al.||2005a^, suggestive of a background 
AGN with intrinsic absorption ( Grindlay et al.|2002{ [Heinke et al.[ 
[2005a^ . Should these identifications as background AGN prove cor¬ 
rect, these sources could be used as calibrators to define a stationary 
reference frame for astrometric studies, allowing an improved dis¬ 
tance determination for the cluster. 

We compared the positions of these three X-ray sources 
(W132, W186 and W231) from the Chandra catalogue of|Heinke et[ 
1^^2005^ with those of their ATCA radio counterparts (Table [^, 
and found that the radio and X-ray frames could be matched to a 3cr 
limit of better than 0"2-0"3. A comparison of our radio position for 
X9 with the X-ray position measured by Chandra in 2002 ( [Heinke[ 
et al.||2005^ suggests that the two are offset by O'.'ll ± O'.'ll in 


R.A. and —0. 12 ± 0. 17 in Dec., i.e. consistent to within Icr (see 

Fig-13'' 


4 47 TUC X9 

An X-ray source was first detected in the core of 47 Tuc by the 
Einstein satellite, with a variable X-ray luminosity that reached 


^ http://www.atnf.csiro.au/research/pulsar/psrcat 
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Figure 3. 5.5-GHz ATCA image of the field surrounding 47 Tuc, encompassing the 8f9 full width at half power (FWHP) of the primary beam at 5.5 GHz (at 
9.0 GHz, the FWHP is 5f3). Red circle shows the 2A" core radius, with the cluster centre marked by the red cross. Sources detected by AEGEAN at 5.5 (9.0) 
GHz are denoted by yellow (cyan) circles, except for X9, which is the one orange circle inside the core radius. 


l.Ox lO^'^ergs"^ (0.5-10keV; Hertz & Grindlay 


1983 


Auriere 


|et al.|19 ^. A faint, variable and very blue optical candidate coun¬ 
terpart, denoted VI, was identified by |Paresce et al.|(T992[ >. ROSAT 
observations resolved the cluster into multiple X-ray sources, of 
which one, X9, was coincident with VI l Hasing er et al.|1994{|Ver-| 
|bunt & Hasinger|1998[ l. [Grindlay et J^ ( |200U used Chandra to con¬ 
firm the association of this source (which they label W42) with V1. 
A far-ultraviolet (UV) observation by |Knigge et al.|(2002| > identi¬ 
fied VI as far-UV bright and variable. 


4.1 Optical/Ultraviolet/X-ray properties 

[Paresce et al.H19^ suggested that the broadband optical through 
ultraviolet spectrum of the source was consistent with an accretion 
disc in a quiescent CV, and on the basis of its high X-ray to optical 


flux ratio, suggested that the system could be an intermediate polar 
(a magnetic CV in which the accretion disc is disrupted inside the 
white dwarf magnetosphere, causing the infalling matter to accrete 
along the magnetic field lines). VI is located blueward of the clus¬ 
ter main sequence in both U against (U—V) and V against (V — I) 
colour-magnitude diagramsJEinonds^etaU200^ and shows sig¬ 
nificant optical flickering jEdmondr'er'n!][2003^ , demonstrating 
that there is a substantial contribution to the optical light from the 
accretion disc. The (non-simultaneous) broadband infrared through 
far-ultraviolet spectral energy distribution (SED) could be quali¬ 
tatively fit by an optically-thick model atmosphere of temperature 
12,000 K and effective radius 0.18f?Q, presumably associated with 
the accretion disc jKnigge et al.|200^ . 

Time-series analysis of the optical counterpart by |Edmonds| 
|et al.| ([?003b|l found a marginally-significant 3.5-hour sinusoidal 
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Table 2. Fitted positions and flux densities of the radio sources detected at > 5.5 it by AEGEAN. Quoted uncertainties are the Itr uncertainties on the elliptical 
Gaussian fit parameters. Identifications show the associated X-ray wavdetect identifications (W-labels) from|Grindlay et al.|p001| and|Heinke et al.|p005a|, 
and the radio source numbers (MA-labels) from |McConnell & Ables|f2000) . Upper limits are given at the 5.5(7 level. The upper limits at 9 GHz vary depending 
on the primary beam response at the source position. 




RA 

(J2000) 



Dec 

(J2000) 


5.5-GHz peak 
flux density 
(^Jy bm~^) 

9.0-GHz peak 
flux density 
(^Jy bm~^) 

ID 

Offset from 

cluster centre 
(arcmin) 

00 

23 

28.608 

± 

0.008 

-72 

03 

53.62 

± 

0.16 

32 

± 

4 

< 76 


3.01 

00 

23 

35.756 

± 

0.007 

-72 

01 

26.24 

± 

0.13 

53 

± 

5 

< 191 

MA2 

4.14 

00 

23 

40.868 

± 

0.002 

-72 

05 

11.62 

± 

0.04 

116 

± 

3 

52 ± 2 

MA3 

1.93 

00 

23 

46.130 

± 

0.008 

-72 

07 

15.21 

± 

0.16 

31 

± 

4 

< 66 


2.81 

00 

23 

52.358 

± 

0.008 

-72 

03 

56.36 

± 

0.15 

27 

± 

3 

< 37 


1.39 

00 

24 

04.247 

± 

0.006 

-72 

04 

58.09 

± 

0.12 

34 

± 

3 

30 ± 2 

W42 (X9) 

0.14 

00 

24 

08.583 

± 

0.009 

-72 

07 

08.40 

± 

0.15 

29 

± 

3 

< 53 

W186 

2.27 

00 

24 

10.304 

± 

0.004 

-72 

07 

41.30 

± 

0.08 

57 

± 

3 

< 66 

MA6 

2.83 

00 

24 

13.313 

± 

0.003 

-72 

03 

32.73 

± 

0.05 

71 

± 

3 

35 ± 1 

MA7 

1.46 

00 

24 

18.081 

± 

0.009 

-72 

06 

19.54 

± 

0.13 

26 

± 

3 

< 44 

W231 

1.73 

00 

24 

27.990 

± 

0.007 

-72 

06 

03.82 

± 

0.14 

30 

± 

3 

< 48 


2.09 

00 

24 

35.034 

± 

0.008 

-72 

05 

01.71 

± 

0.15 

33 

± 

4 

< 48 

W132 

2.26 

00 

24 

47.350 

± 

0.010 

-72 

07 

53.63 

± 

0.12 

40 

± 

4 

< 249 


4.40 

00 

24 

55.057 

± 

0.010 

-72 

05 

28.50 

± 

0.14 

37 

± 

4 

< 144 


3.84 


Table 3. Fitted 5.5-GHz positions and flux densities of the MSPs detected at > 3cr. All known MSP positions were inspected by eye. Where radio emission 
was detected at the pulsar position, a point source fit was caiTied out using the AIPS task JMFIT. 


RA 

(J2000) 

Dec 

(J2000) 

Peak flux density 
(/iJy bm~^) 

ID 

00 23 59.49 ± 0.04 

-72 03 58.8 ± 0.2 

19 ±5 

MSP-J 

00 24 04.68 ± 0.05 

-72 04 54.4 ± 0.3 

16 ±5 

MSP-0 

00 24 06.71 ± 0.07 

-72 04 06.7 ± 0.4 

11 ± 5 

MSP-H 

00 24 08.02 it 0.05 

-72 04 39.6 ± 0.3 

15 ±5 

MSP-G/I 

00 24 08.47 ± 0.05 

-72 04 38.6 ± 0.3 

14 ± 5 

MSP-T 

00 24 13.92 ±0.03 

-72 04 43.4 ± 0.2 

23 ±5 

MSP-D 

00 24 16.63 ± 0.05 

-72 04 25.0 ± 0.3 

16 ±5 

MSP-Q 


signal in the U-band power spectrum that implied an orbital pe¬ 
riod of ~ 3.5 or 7 hr, depending on whether or not the periodic¬ 
ity could be ascribed to ellipsoidal modulations. For a Roche-lobe 
filling main-sequence donor star, |Knigge et al.| ( |2008[ l found that 
a 7-hr period would imply a donor that was too bright at wave¬ 
lengths > 5000 A. For non-ellipsoidal variability, a 3.5-hr period 
would then imply a donor of mass 0.5Mq, which, while not re¬ 
quired, could be accommodated by their model for the broadband 
spectrum, which was dominated by the hot accretion disc. In slit¬ 
less ultraviolet spectroscopy with HST, |Knigge et al.| ( |2008[ l also 
detected the spectroscopic signature of a relatively high-inclination 
accretion disc or disc wind, namely a strong, double-peaked C IV 
emission line, which they took to confirm the identification of the 
source as a CV. 

With a 0.5-6 keV luminosity of 0.6-l.Ox 10^® erg s“^ 

|et al.|2005a[ after scaling to our adopted distance of 4.57 kpc), the 
high X-ray-to-optical flux ratio (Fx /Fopt ~ 5) of VI is consistent 
with the system being a CV, but is higher than any known field CV, 
and comparable only to another (confirmed magnetic) CV in 47 


Tuc, V3/X10 I 

Edmonds et aL|20U3b(|Heinke et al.|2005a|l. While 

[Grindlay et al. 

pool [I reported a marginal detection of a 218-s pe- 


riodicity in the X-ray emission from X9, the original claim was not 
highly significant, and attempts with larger datasets to recover this 
periodicity have failed (Bahramian et al. in prep.). 


The 2000 Chandra X-ray spectrum of X9 could be fit by a 


Heinke 


thermal bremsstrahlung model with a temperature of 32 ± 6keV 
( [Grindlay et al.| |2001^. However, a higher-quality 2002 Chandra 
spectrum revealed additional features ( [Heinke et al.|[2005a| l. A 
strong emission line (likely O VIII) at 0.65 keV superimposed on 
a hard continuum could be fit with a thermal plasma model, with 
two MEKAL components at 0.25 and > 17keV, together with a 
third, highly-absorbed component needed to fit the spectmm above 
6 keV. This is not typical of a CV, but equally is not obviously in¬ 
dicative of any other class of source. However, the long exposure 
time (~ 70 ks in 2000 and ~ 280 ks in 2002) and the low ab¬ 
sorption towards 47 Tuc (A^h = 1-3 x 10^°cm“^; 

|2005a^ imply that the spectrum is of higher quality than typically 
available for accreting stellar-mass compact objects. Since fits to 
MEKAL components with high temperatures are typically indistin¬ 
guishable from power-law spectral fits, we therefore re-fit a short, 
5.3-ks observation (ObsID 3384, from 2002 September 30, taken in 
subarray mode to alleviate pileup) with a model consisting of an ab¬ 
sorbed power-law plus a Gaussian line to represent the low-energy 
MEKAL component, to enable a comparison with typical spectra 
from other classes of source (see Sections [4.3| and [4!^ . Fixing the 
absorption to the cluster value, we found a fitted photon index of 
F = 1.06±0.24, aline energy of 0.56±0.03keV, and aline width 
of0.10±0.03keV. 


Heinke et al. 
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4.2 Radio emission from cataclysmic variables 


Our radio detection of X9 provides new constraints on the nature 
of the accreting object. Since the source was previously identified 
as a CV, we first compare its radio properties with those of known 
CVs, of which only a handful have been detected in the radio band 
to date. 


Of the magnetic CVs, the most luminous persistent radio emit¬ 
ter is the rapidly-rotating intermediate polar AE Aqr, with a typical 
5-GHz flux density of 4-8mJybeam“^. Its spin period of just 33 s 
I Patterson|1979) creates a p ropeller effect jWynn, King & Home| 
r99^| ~ |Eradeous & Horne||1996^ , with the expulsion of magne¬ 
tised plasma clouds giving rise to the observed, rapidly-variable 
synchrotron emission jMeintjes & Venter|2005^ , which has a typ¬ 
ical spectral index of 0.35-0.60. At a distance of 86 pc, its max¬ 
imum measured 5-GHz flux density of 16mJy (during a short- 
duration, ~ 2-hour flare; see|Bastian, Dulk & Chanmugam| 1988[ 
|Abada-Simon et al.|1993^ would correspond to a radio luminosity 
of 7 X 10^“ergs“^ ' 

Other than AE Aqr, only the two polars AM Her (|Chanmugam| 
|& Dulk|1982| l and AR UMa ( [Mason & Gray|200'7) have been found 
to be persistent radio emitters among the magnetic CVs. Their qui¬ 
escent emission (Lr ~ 2 x 10^® erg s“^) has been attributed to gy- 
rosynchrotron emission from energetic electrons in the white dwarf 
magnetosphere, with occasional highly circularly-polarised flaring 
events due to cyclotron maser emission ( [Mason & Gray||2007[ (. 
The quiescent emission of AM Her appeared to be relatively sta¬ 
ble over a ten-year period jChanmugam & Dulk[|1982| [Ravelin, [ 
[Spencer & Davis|1994^ , and from observations at 4.9 and 8.4 GHz 
taken 22 years apart, we derive a non-simultaneous spectral index 
of-0.3 ±0.3. 


The most luminous reported radio emission from a CV was 
a decaying flare detected from the intermediate polar DQ Her 
by [Ravelin et al.[jl994^ . Over 30 min, the flux density decreased 
by a factor of 4, with an initial radio luminosity of = 5 x 
10^^ergs“^. The system was not detected on three further occa¬ 
sions. Similarly short-timescale flares with comparable luminosi¬ 
ties have also been reported from V834 Cen jWright et al.|198^ 
and BG CMi jRavelin et al.|1994l l, although these detections were 
both lower-significance, and the former was from a single-dish ob¬ 
servation and therefore potentially more susceptible to the effects 
of both source confusion and radio frequency interference. 

Of the non-magnetic CVs, the only published radio detections 
to date are from the dwarf novae SS Cyg, EM Cyg and SU UMa, 
and the nova-like systems V3885 Sgr, TT Ari, RW Sex and V603 
Aql. Of these, the most well-studied is SS Cyg, whose radio coun¬ 
terpart has been repeatedly detected at the peak of its frequent out¬ 
bursts, and whose high brightness temperature radio emission has 
been interpreted as synchrotron emission from a jet ( [Kdrding et al.[ 
[2008[ l. At a distanc e of 114 ± 2 pc, its brighte st detected radio flare 
of 3mJybeam~^ jMiller-Jones et ah 2013 ' corresponds to a ra¬ 
dio luminosity of 2.3 x 10^®ergs“^. Single detections of both SU 
UMa ( [Benz et al.[1983^ and EM Cyg ( [Benz & Gudel[198^ implied 
similar radio luminosities of ~ 1-5 x 10^®ergs“^ (depending on 
the poorly-constrained source distances). The four no va-like sys- 
tems were detected at levels of 0.03-0.2 mjybeam“^ jKording'rt 


[al.|201 l1[Coppejans et al.|2015[ l, which for the estimated distance 
ranges imply radio luminosities of 4.5 x 10^'*-!.6 x 10^® erg s“^. 


In summary, our radio detection of X9 is an order of magni¬ 
tude more luminous than the persistent emission detected from any 
known CV. While it is comparable to the peak luminosities of in¬ 
frequent short-timescale flaring events, the stability of the source 


over the course of each 10-h observing run, and between the two 
epochs in 2010 and 2013, argues against an origin as flaring events. 
Furthermore, the absence of a high level of circular polarisation 
suggests that the emission is unlikely to arise from an electron cy¬ 
clotron maser. Thus, in light of our radio data, we do not favour an 
interpretation of X9 as a C V. 


4.3 A transitional millisecond pulsar? 

Another class of accreting sources known to produce flat-spectrum 
radio emission are the binary MSRs that show occasional transi¬ 
tions to and from an accretion-powered state characteristic of an 
LMXB. To date, we know of three such confirmed transitional bi¬ 
nary systems, which can spend years in a phase of quiescent ac¬ 
cretion, with X-ray luminosities in the range erg s“^ ( [d^ 

[Martino et al.[[2013| [Ratruno et al.[[2014) . Both in this quiescent 
phase and also during the brighter outburst phase the observed mil¬ 
lisecond X-ray pulsations at the pulsar spin period show that some 
fraction of the accreted material is reaching the neutron star sur¬ 
face ( [Rapitto et al.[2013][Archibald et al.[2015[[Rapitto et al.[2015^ . 
All three systems have also shown flat-spectrum continuum radio 
emission when in their accretion-powered states ^Bond et al.[2002[ 


Rapitto et al.|2013^[Hill et al.|201 l[|Bassa et al.|2014[ Deller et al.[ 


2015 I, with inferred radio luminosities significantly higher than ex¬ 
pected from an extrapolation of the radio/X-ray correlation for the 
traditional atoll source neutron stars (Fig.|^ see [Deller et al.[20T5) l. 
However, the proposed correlation between their measured radio 
and X-ray luminosities ( [Deller et al.|2015) implies radio emission 
a factor of ~ 6 lower than seen for X9 at a similar X-ray luminosity 
(see Fig.|^. 

The three transitional MSR systems share some unusual X- 
ray variability characteristics, which can be compared to the X- 
ray variations of X9. Each transitional system has switched at 
least once between a radio pulsar state (with Lx ~ 10®^ erg s“^) 
and a low-luminosity accretion state (with Lx averaging ~ 4 x 
10®® erg s“®), with M28I also showing a bright accretion outburst 
at Lx > 10®® ergs“®. Within the low-luminosity accretion state, 
each system switches among three modes; a “high” mode showing 
pulsations, sudden dips to a “low” mode (transitions less than 30 s, 
average flux a factor 5-10 lower), and a rarer “flaring” mode (last¬ 
ing 1-45 minutes, typically a factor >4 brighter). This state also 
shows a remarkably constant average X-ray flux (set by the “high” 
mode flux); multiple observations over years have not found factor- 
of-two variations jde Martino et al.[2013[[Linares et al.[2014[[Bog-[ 
[danov et al.[2015[ l. This is in contrast to X9, where the average flux 
varies smoothly by a factor of ten over timescales of several days 
(Bahramian et al. 2015, in prep.). 

Each transitional MSB system also shows all three modes. 
Dips to the “low” mode can be quite short (< 100 s), below our ca¬ 
pacity to detect in the lightcurve of X9. However, dips longer than 
100 s occur frequently; in RSR J1023-I-0038, dips longer than 200 s 
occur every 3 ks ( [Bogdanov et al.[2015) , and in XSS J12270—4859, 
dips between 100 and 1100s also occur on average every 3ks 
[Martino et al.|2013y The M28I lightcurve of [Linares et al.[ ( [2014| , 
limited to binning of 500 s, does not probe such short timescales, 
but nevertheless showed six low states of length exceeding 1 ks, 
spread over 234 ks. Flares are also common, occurring every 20 ks 
or so in RSR J1023-I-0038 ( [Bogdanov et al.[2015[ l, and at similar fre¬ 
quency (occupying 1/10 to 1/3 of the time) in XSS J12270—4859 
([de Martino et al.[2013) . Careful inspection of the 300 ks of ACIS-S 
lightcurves of X9 from 2002 ( [Heinke et al.|2005^ , binned at 50 s, 
fails to show evidence of either sharp transitions to low modes, or 
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Figure 4. Radio/X-ray correlation in hard state and quiescent stellar-mass black holes. Black squares show field black holes from the literature jMiller-Jones 


et al. 

2011 

et al. 

2011 

Gray 

2007 


Papitto et al.|2013||De ler et al.|2015), green diamonds show the magnetic CVs jAbada-Simon et al.|1993||Chanmugam & Dulk|1982l|Mason & 

yellow circle shows the M22 candidates jStrader et al.|2012a) , red circle shows the M62 candidate' 7Qiomiuketari2013) , and purple circle shows 
47 Tuc X9. Filled symbols are from simultaneous radio and X-ray data, where as open symbols are non-simultaneous measurements. Black dotted line shows 
the best fit relation for the black hole systems (Lr oc Gallo et al. 2006 . Blue dashed and dash-dotted lines show two possible radio/X-ray correlations 


for neutr on star system s ^Migliari & F ender|200^ , and green line shows the maximum radio luminosity of non-magnetic accreting white dwarf systems in 
outburst jKdrding et al.|2008]|20TT^. X9 falls within the scatter of the best fitting radio/X-ray correlation for black holes, and is too radio luminous to be a CV. 


flaring by a factor >2. Extrapolation from the transitional MSPs 
suggests that of order 100 low modes, and 15 flaring episodes, 
should have been observed; the lack of these behaviours strongly 
suggests that X9 is not a similar system as the transitional MSPs. 

Furthermore, the optical spectrum of PSR J1023-I-0038 dur¬ 
ing its accretion-powered state shows multiple double-peaked He 
emission lines ( |Wang et al.|2009^ , in contrast to X9, which shows 
no strong evidence for He emission in its far-ultraviolet spectrum 
( [Knigge et al.|2008^ . The X-ray spectrum of X9 (Section [4T| l is also 
markedly different from the spectra of the transitional MSPs when 
at comparable luminosities in the accretion-powered state. With a 
similarly low absorption column density of Wh ~ 3 x cm“^, 
PSR J1023-I-0038 shows no evidence for an emission line near 
0.6 keV, and the spectrum can be well fit by a simple power-law 
model with no requirement for an additional MEKAL component 
( [Bogdanov et al.|2015[ l. And despite the larger column densities to¬ 
wards XSS J12270—4859 and M281, there is no evidence for such 
an emission line in high-quality spectra of either source jde Mar-j 
jtino et al.|2010[ [Papitto et al.|2013| l. Furthermore, at X-ray lumi- 
nosities of erg s“^ all three transitional MSPs can be fit 

by simple power-law spectra with photon indices of F =1.5-1.7 
( [de Martino et al.|2010[[Papitto et al.[2013[[Bogdanov et al.[2015[ l, 
significantly softer than seen in X9 (Section fT.l^ . 

In summary, while we cannot explicitly rule out X9 being a 


transitional MSP, the available X-ray, optical and ultraviolet data, 
taken together with the relatively high radio luminosity, lead us to 
consider this explanation to be unlikely. 

4.4 X9 as a black hole 

The high radio luminosity implied by our observations of X9 sug¬ 
gests an alternative identification of X9, as a stellar-mass BH. Com¬ 
parison of our radio measurements with the X-ray luminosity mea¬ 
sured by [Heinke et al.[p005a[ l shows that the source lies very close 
to the well-established radio/X-ray correlation for accreting stellar- 
mass BHs (Fig. 1^. And while the X-ray and radio observations 
were separated by 8-10 years, the relative stability of the radio flux 
density between 2010 and 2013, and the X-ray flux between 1992 
and 2002 jVerbunt & Hasinger|1998[[Grindlay et al.|2001[[Heinke] 
[et al.[2005a| l helps assuage (but does not eliminate) concerns about 

strong source variability. 

[Knigge et al.[ p008[ l found the broadband optical-through- 
ultraviolet spectrum of V1 to be dominated by an optically-thick 
component of temperature ~ 12,000 K and radius ~ O.ISJJq, con¬ 
sistent with a disk-dominated nova-like CV. However, as noted 
by [Marsh, Robinson & Wood| ( [1994l l, there are no strong differ¬ 
ences between the optical spectra of dwarf novae and those of qui¬ 
escent black holes. In the case of VI, we note the similarity to 
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the broadband spectra of several dynamically-confirmed quiescent 
BHs I McClintock et al.|2003[[Froning et al.|201 l[[Hynes & Robin-| 

|son l|2012t , which could all be fit by blackbodies of temperature 

5,000-13,000 K, with areas much smaller (0.06-1.25%) than the 
projected area of the expected truncated quiescent accretion disk 
(i.e. radii of 0.01-0.09J?q). This emission was explained as either 
the impact point of the accretion stream on the disc ( [McC lintock et| 
|al.|2003f, or as emission from the inner part of the disc, close to the 
truncation radius ( [Hynes & Robinson|2012^ . 

The main caveat to this explanation is that the X-ray spectrum 
of X9 is difficult to reconcile with those of known quiescent black 
holes, which at luminosities of 10®^ erg s“^ are fit by relatively soft 
power-law spectra with photon indices of T ~ 2.1 [Plotkin et al.[ 
[2013[ l. While no evidence for low-energy emission lines close to 
0.6 keV has been observed in any quiescent black hole, the high col¬ 
umn densities towards most such systems would have precluded the 
detection of such a line. Only three known quiescent systems (XTE 
J1118-H480, A0620-00 and MAXI J1659-152) have a bsorption 
column densities < 5 x 10^^ cm“^ ( Plotkin et al.|2013) , none of 
which have stringent published constraints on the presence of low- 
energy emission lines. The only published quiescent spectrum of 
sufficiently high quality to investigate the presence of such lines is 
that of V404 Cyg ( [Bradley et al.[2007l l. That spectrum was fit with 
an absorption column of 7-9 x 10^^ cm and clearly excludes a 
line of similar strength in the range 0.5-0.7 ke\|^ In summary, the 
X-ray spectrum, while not fully characteristic of a quiescent black 
hole, is not better explained by any other known class of source. 


5 THE NATURE OE THE DONOR 

Having identified a black hole as a plausible accretor, we now ex¬ 
plore the properties of X9 to help determine the nature of the sys¬ 
tem. 


5.1 Orbital period 


We reanalysed the optical data of [Edmonds et al. [ ( [2003b) to check 
the significance of the marginal 3.5-hour sinusoidal signal and 
found the false detection probability to be > 10 per cent, imply¬ 
ing a significance of < 2a. We also searched for shorter periods 
than the 1.5-14hour range considered by [Edmonds et al.][2003b[ l, 
finding nothing significant at the 3a level. Thus we have no good 
photometric constraint on the orbital period. 

The measured U-magnitude of X9 is in the range 19.8-21.4 
( [Edmonds et al.|2003a[[Beccari et al.|2014^ , which at 4.57 kpc im¬ 
plies an absolute magnitude in the range 6.5-8.1. Should the op¬ 
tical light be attributed to the donor star, this rules out a subgiant 
or high-mass companion and implies a lower main-sequence donor 
of spectral class later than mid-K. Should the optical emission be 
dominated by the accretion flow (as seems likely given the spec¬ 
tral fitting of [Knigge et alj 2d0^ the donor star would be fainter 
still. As determined by Knigge et al.[ ( [2008) l, periods shorter than 
7 hr would be required for Roche-lobe overflow of a main-sequence 
donor star with a sufficiently faint visual magnitude. 

At orbital periods longer than a few hours, mass transfer in 


^ The published spectrum shows residuals close to 0.8 keV, but they are not 
statistically significant, and a line at this energy has no clear explanation, in 
contrast to the 0.65 keV line in X9 that could potentially be explained as 
O VIII. 


X-ray binaries is driven by magnetic braking or nuclear evolution 
of the donor star. In this regime, there is a well-established posi¬ 
tive empirical correlation between orbital period and quiescent X- 
ray luminosity ( Gareia^el^^y[2TOl|| McClintock, Narayan & Ry-| 
[bicki[[2004{ [Arrnas Padilla et al.[[2014]f 7 as expected from binary 
evolution theory in the case where quiescent accretion proceeds 
via an advection dominated accretion flow ( [Menou et al.[[T99^ . 
At shorter orbital periods (less than a few hours), angular mo¬ 
mentum loss via gravitational wave radiation will take over as the 
dominant mechanism driving the mass transfer. This should lead 
to mass-transfer rates (and hence X-ray luminosities) that increase 
with decreasing orbital period, implying the existence of a mini¬ 
mum X-ray luminosity at a bifurcation period of ~ 10 h ( [Menou eT| 
[al.|1999[ l. While there is a dearth of known black hole systems at 
short orbital periods, there is some evidence for a minimum qui¬ 
escent luminosity [Gallo et al.[[2008% at a level consistent with 
the best current constraints on the recently-discovered BH candi¬ 
dates in the GC M22 ( [Strader et al.[2012^ . The recent detection of 
two black holes with orbital periods < 3 h has begun to probe the 
short-period regime, with MAXI J1659—152 Homan et al.|2013[ l 
and Swift J1357.2—0933 ( [Armas Padilla et aE 2014[ l both showing 
higher quiescent luminosities than systems with orbital periods a 
few times longer (although the uncertain distance prevents a defini¬ 
tive conclusion in the latter case). 

At 10^^ erg s“^, the X-ray luminosity of X9 is high for a qui¬ 
escent BH, exceeding that of V404 Cygni (Eig. [^, the most lu¬ 
minous known quiescent field BH X-ray binary. The empirical re¬ 
lationship between orbital period and quiescent X-ray luminosity 
would then imply either a period that is roughly comparable to the 
6.5-day period of V404 Cyg ( [Casares et al.|1992l l, or else a very 
short orbital period in an ultracompact system (i.e. with a white 
dwarf or helium star donor). The former scenario is ruled out by the 
measured optical emission, although given the scatter on the empir¬ 
ical period-luminosity correlation, we cannot exclude the possibil¬ 
ity of a low-mass main-sequence donor with a period < 7 h. Should 
the system be ultracompact, the absence of He II 1640 A line emis¬ 
sion in the far-UV spectrum of [Knigge et al.[ ( [2008[ would suggest 
a white dwarf rather than a helium star donor. 

Although no Galactic ultracompact BH systems have yet been 
discovered, it is possible to form such systems in a GC environ¬ 
ment following an exchange interaction, a long phase of subsequent 
hardening through stellar encounters, and finally triple-induced 
mass transfer ( [Ivanova et al.[[20T0) . Indeed, the black hole in the 
extragalactic GC RZ2109 in NGC 4472 is believed to be an ultra¬ 
compact system with a white dwarf donor and an orbital period of 
~ 5 min ( [Zepf et al.[2007[|2008[[Gnedin et al.[2009[ l and also shows 
an excess of soft photons that can be modelled as O VIII emission 
( [Joseph et al.|2015) l, so such a scenario would not be wholly un¬ 
precedented. 

Given the similarities between the quiescent luminosities of 
X9 and V404 Cyg, w e might exp ect similar mass transfer rates, 
i.e. ~ 10“®Mq yr-i |King|l9^ . Such a high mass transfer rate 
would require the system to be ultracompact, with an orbital pe¬ 
riod of ~ 20 min ( [van Haaften et al.|2012[ . While neutron stars in 
ultracompact systems with such short orbital periods are all persis¬ 
tent X-ray sources ( [Heinke et al.|201^ , the low X-ray luminosity 
of X9 implies that it must be a quiescent transient system, which 
undergoes occasional outbursts due to a disc instability. Since the 
critical mass transfer rate for the disc instability model scales with 
accretor mass as then black hole systems will become tran¬ 

sient at slightly higher mass transfer rates. [van Haaften et al.[p0T^ 
see their fig. 8) show that the highest mass transfer rate for which 
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a black hole ultracompact system will become transient (a few 
tens of Myr following the onset of mass transfer) is a few times 
yr“^ (assuming a pure He donor), with an orbital period 
of ~ 25 min and a donor of mass 0.02Mq . If the mass transfer rate 
is to remain approximately similar to that of V404 Cyg, the system 
should not be significantly older (and with a consequently longer 
orbital period) than this. 

For the disc instability model to operate, the outer disc would 
need to be below the ionisation temperature of the material, imply¬ 
ing that the ~ 0.187?q optically-thick component fit to the SED by 
|Knigge et al.| ( |2008) l can not extend over the entire disc. Should the 
system he ultracompact, then for a 25-minute orbit and a BH mass 
of 3-20Mq, the circularization radius should be 1.5-3.7 x 10^° cm 
(from lron Haaften et al.|2012| fig. 5), so the fitted emission size 
would be 34-85% of the circularization radius. While this might 
seem uncomfortably large, the physical size is comparable to (albeit 
slightly larger than) that seen in other quiescent BH X-ray binaries 
([McClintock et al.|2003[ [Froning et al.|2011[|Hynes & Robinson| 


|2012 as discussed in Section |4.4^ , and is consistent with typical 

truncation radii of quiescent discs (a few thousand Schwarzschild 
radii; see, e.g. |Menou et al.|2000| l. While a cooler outer disc might 
be easier to accommodate in a longer-period, main-sequence sys¬ 
tem, an ultracompact system remains plausible, particularly given 
that the non-simultaneity of the data and the time-varying nature 
of the source permitted only a qualitative fit to the SED. Quasi- 
simultaneous, broad-band photometry would be required to prop¬ 
erly characterise the size of the emitting region and hence ascertain 
the viability of this scenario. 

Within the disc truncation radius, quiescent black holes host 
radiatively inefficient accretion flows (e.g. advection dominated ac¬ 
cretion flows; ADAFs, [Narayan & Yi|1994) . with a population of 
hot electrons producing the observed X-ray emission via thermal 
Comptonisation. In an ultracompact system, the accreted material 
would consist of more massive ions than the hydrogen-rich material 
being accreted by all known quiescent black holes in the field. The 
additional gravitational energy liberated could potentially increase 
the temperature of the particles within the ADAF, thereby affecting 
the photon index of the power-law and leading to a harder spectrum 
for X9 in comparison to known quiescent black holes. However, an 
in-depth exploration of this issue is beyond the scope of this paper. 

5.2 Constraints on the Ha emission 

Since an ultracompact system must have a hydrogen-deficient, de¬ 
generate donor star, a test for the presence of Ha emission could 
also help to confirm or refute the ultracompact scenario. Whereas 
no Balmer emission would be expected from an ultracompact 
system, the Ha equiva l ent wi dth (EW) of V404 Cyg was mea- 


analyzed an HST data set of quasi-simultaneous images of 47 Tuc 
in Ha and the surrounding continuum. The observations were taken 
with the Wide Field Channel (WFC) on the Advanced Camera 
for Surveys (ACS) under program GO-9281 (PI: J. Grindlay). Im¬ 
ages were obtained through blue (F435W), red (F625W) and Ha 
(F658N) filters, with total exposure times of 955 s, 1320 s, and 
7440 s, respectively. The observations were split over three vis¬ 
its, on 2002 September 30, October 2, and October 11, with ob¬ 
servations being taken in all three filters at each epoch. The data 
were taken simultaneously with the 2002 Chandra observations of 
47 Tuc presented in |Heinke et al.|p005a| l; a full description of the 
GO-9281 data and results is forthcoming (van den Berg et al. in 
preparation). 

We retrieved the flat-fielded images, corrected for charge- 
transfer efficiency losses, from the HST archive. These images still 
suffer from significant geometric distortion, which we corrected for 
with the DrizzlePac software jGonzaga et al.|20T^ . The individual 
images were astrometrically aligned and subsequently combined 
for each filter using astrodrizzle, which creates stacked images that 
are cleaned of cosmic rays, and that are without geometric dis¬ 
tortion. These master images are twice-oversampled to 0^025 per 
pixel. All images except two 10-s F625W exposures and one 10-s 
F435W exposure, were included in the stacks. 

Photometry for X9 and the surrounding stars was extracted us¬ 
ing the point-spread-function fitting package DAOPHOT. The mag¬ 
nitudes were calibrated to the VEGAMAG system using the zero- 
points on the STScI websit^ The position of X9 in the colour- 
magnitude diagram (Fig.|^ shows that it is indeed a very blue ob¬ 
ject, but, in contrast to the findings of |Beccari et al.| ( (2014| l, it does 
not stand out in (Ha — R% 25 ) colour from other stars of similar 
7?625 magnitude. Since the average (Ha — Rs 2 s,) colour of stars 
varies with (H 435 — Rq 25 ) colour, we also compared the position 
of X9 in a colour-colour diagram (Fig.|^. Unfortunately, X9 is one 
of the bluest stars in the cluster, so there is a dearth of similarly blue 
stars. Instead, an estimate of the Ha equivalent width can be made 
by comparing the unreddened colours of X9 with those expected 
for the blue component believed to dominate the spectrum of X9 
between ultraviolet and infrared wavelengths ( [Knigge et al.|2008) ). 

If X9 is an ultracompact binary, a 12,000 K blackbody would 
be a reasonable approximation of the optical spectrum. The differ¬ 
ence between the actual and expected (Ha — Rs 2 z>) colour then 
suggests an excess emission in the Ha filter that is similar to the 
effect of adding an Ha emission line with an EW of about —10 A. 
In an ultracompact binary, no Ha emission is expected. The small 
inferred excess could point at an incorrect description of X9’s con¬ 
tinuum emission by the assumed blackbody, or perhaps at the pres¬ 
ence of an emission line other than Ha; for example, the CII/O II 
emission line at 6580 A in the spectrum of t he ultracompact binary 


sured by jCasares et al. 1 1993 1 to be —38.7 A. This is roughly as 4U 0614-1-091, which has anEW of ~ —3 A 1 Nelemans et al. 


2004 


expected from the anticorrelation between X-ray luminosity and 
EW reported by [Fender et al.| (|2009[), for a system of luminosity 


[Baglio et al.[[2014 

1 . Since the X-ray luminosity of 4U 0614-1-091 

(3 X 10^“ergs“^ 

Migliari et al. 12010 1 is significantly higher than 


10 ^® erg s"^ 


[Beccari et al.j \2QIA) used archival HST data to calculate 
(Veoe— Ha) and (V606-7si4) colours for X9, which showed an Ha 
excess with an equivalent width of —125 A. However, their Veoe 
and 7 gi 4 colours on the one hand, and Ha data on the other, were 
taken in different years, and a comparison of their V -band (F606W) 
data from 2006 March with F555W observations from 1999 July 
( [Edmonds et al.|2003a^ shows a difference of 1.55 magnitudes, suf¬ 
ficient to render their non-simultaneous (Veoe— Ha) colour mean¬ 
ingless. 

To obtain a better measurement of the Ha emission in X9, we 


that of X9, the outer disc in X9 is expected to be more optically 
thin in the continuum. Since the line emission remains optically 
thick (e.g. |Williams|1980l (, the EW of these C/O lines might be ex¬ 
pected to be larger in X9 (as also manifested in the anticorrelation 
between X-ray luminosity and Ha EW; [Eender et al.|2009^ . 

If X9 is not an ultracompact binary, then the appropriate com¬ 
parison model should include absorption lines as seen in the spec¬ 
tra of optically thick accretion disks. Comparing the colours of X9 


http://www.stsci.edu/hst/acs/analysis/zeropoints 
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435“ ^625 Ha-Rg25 

Figure 5. Colour-magnitude diagrams showing the -B 435 , -R 625 , and Ho photometry obtained from the HST ACSAVFC GO-9281 images of 47 Tuc. The 
optical counterpart to X9 is shown in red, and the surrounding stars are plotted with black points. For comparison, we also show (in red) the positions of the 
optical counterparts to three other Chandra sources in 47 Tuc: V3 (W27/X10) and W45 (classified as CVs) and X5 (W58; classified as a quiescent NS LMXB). 


with those for an early A-type star would then suggest an Ha EW 
of about —20 A. In this case, part of the excess emission could be 
explained by a true Ha emission line associated with the accretion 
disk, or part could be emission from the donor star, or perhaps a 
disc wind. The observed flux at red wavelengths in the SED of X9 
could accommodate the signature of a cool (~ 4000 K) secondary 
( [Knigge et al.|2008l i, but not that of an A-type main sequence star. 
Furthermore, given that X9 is variable and that the data points that 
make up the SED were not taken simultaneously, such an additional 
component may not be required. Regardless, it is clear that there is 
no strong evidence for an Ha excess from X9. 

We can also compare our estimate of the Ha EW in X9 
with that found by [Wang et aT] ( |2009| l for the transitional MSP 
PSR J1023-b0038. Summing both the red and blue components of 
the Ha line in PSR J1023-I-0038, |Wang et al.| found a line flux 
of 1.8 X 10“^'* ergcm“^ s“^. Combining the i ?625 magnitude of 
19.8 with our estimated Ha EW for X9 of —10 A and account¬ 
ing for the difference in source distance (assuming 1.37 kpc for 
PSR J1023-I-0038; [Deller et al.|[M^ , we find the Ha line lu¬ 
minosity of X9 to be a factor of ~ 7 times lower than that of 
J1023-I-0038, suggesting that X9 has either a smaller disc or less 
Ha emission. Despite similar absolute magnitudes in the U band 


HKnigge et al.||2008| [PrEuno et al.|[M4t , PSR J1023-I-0038 has 
a redder continuum than X9 (a {B — R) colour of 1.0 compared 
to 0.4 for X9). A comparison of the optical magnitude of the ac¬ 
creting (disc-dominated) state of PSR J1023-I-0038 with that in its 
rotation-powered state suggests that the optical continuum emis¬ 
sion is dominated by the disc rather than the donor star in the 
accretion-powered state (| Wang et al.||2009l [Thorstensen & Arm-| 
|strong|2005| l. Thus, the redder continuum implies more emission 
from the outer disc. Together with the correlation between orbital 
period and quiescent X-ray luminosity, the anticorrelation between 
Ha EW and quiescent X-ray luminosity ( [Fender et al.|200^ im¬ 
plies that shorter-period systems (with a smaller disc) should have 
a higher Ha EW. Since the larger disc in PSR J1023-I-0038 has a 
higher Ha EW, we therefore conclude that X9 is likely hydrogen- 
deficient. 

Finally, we note that ~ 70% of CVs are found to be Ha emit¬ 
ters, with the only magnetic systems found to lack Ha emission 
being polars in the low state jWitham et al.|2006l l. This provides ad¬ 
ditional circumstantial evidence against the CV scenario, although 
we note that the sample of jWitham et al.j included only a single 
intermediate polar. 
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Figure 6. Colour-colour diagram showing in red the optical counterparts 
to X9 and the three comparison Chandra sources from Fig. Colours 
of suiTounding 47 Tuc stars are shown as black points. The photometry 
has been dereddened to account for the small extinction towards 47 Tuc 
(E{B — V) = 0.04). The expected colours for main-sequence stars of 
various spectral types are shown in blue, and were computed with the 
synphot package using low-metallicity (log Z = —2.5) Castelli-Kurucz 
model spectra. The solid blue line is for EW(Ho)=0 A, and the dotted lines 
show the effect of adding a gaussian emission line with a FWHM of 25 A in 
increments of AEW(H«)=—10 A. In green we show the synthetic colours 
for a 12,000 K blackbody model jKnigge et al.|2008} with an added gaus¬ 
sian emission line of different EWs. 


6 DISCUSSION 


we determined the number of transient systems of different types: 
29 LMXBs with dynamically-confirmed BHs or strong BH candi¬ 
dates; 14 LMXBs that are weak BH candidates (and thus could be 
confused with NS LMXBs); 60 LMXBs with either Type I X-ray 
bursts or no indication of being a BH (and thus are assumed to be 
NS LMXBs). Ignoring the harder-to-differentiate weak BH candi¬ 
dates, we estimate the fraction of transient LMXBs in the field that 
are BHs to be We can then estimate the probability 

of observing zero BH LMXBs in GCs given the ten known tran¬ 
sient NS LMXBs, assuming the field and GCs share the probability 
density function of a LMXB containing a BH in the field. We find 
this null hypothesis probability to be just 2.3%. While this suggests 
there may be a difference between the field and GC populations of 
transient LMXBs, being < 3 (t it cannot be taken as conclusive ev¬ 
idence. We note that this statement remains true even if we classify 
all the weak BH candidates in the field as BHs. 

Furthermore, we note that three of the four cluster BH candi¬ 
dates identified to date could be ultracompact systems (see |Strader| 
|et al.|2012a| for M22, and Section|^for 47 Tuc), whereas no such 
systems are known in the field. |Knevitt et al.| ( |2014| l suggest that 
BHs in short orbital period binary systems could remain undetected 
due to selection effects. The correlation between orbital period and 
peak outburst luminosity ( |Wu et al.|2010| l implies that the luminos¬ 
ity of short orbital period BHs may never exceed the threshold for 
radiatively efficient accretion, causing a system to remain in a low- 
luminosity hard state throughout its outbursts ( |Meyer-Hofmeister| 
|2004| >. The weak irradiation also leads to significantly shorter out¬ 
bursts, and the duty cycles are similarly reduced to low levels. The 
combination of these effects then implies that such outbursts could 
have been undetectable or overlooked by current or previous All 
Sky Monitors. A more sensitive future wide-field monitor such as 
proposed for LOFT ( [Maccarone et al.||2015| ) would be needed to 
detect any such faint outbursts. 


X9 is the fourth BH candidate in a Galactic GC discovered via deep 
radio continuum observations. By design, the clusters selected for 
our deep radio observations (M22, M62 and 47 Tuc) were all rel 
atively massive (M > 3 x 10® Mq; Pryor & Meylan|l9'^ 


Kim- 


|mig et al.|2015^ . In the absence of significant natal kicks, the large 
masses imply a significant initial population of BHs, and the long 
evolutionary timescales make them more likely to retain BHs jHeg-| 
Igie & Giersz||2014l>. While the absence of any BH candidates in 


M15 I Strader et al.|2012b) might suggest that other factors are at 
play, we note that at 10.4 kpc ( |Harris|1996| l, M15 is over 50% more 
distant than the most distant of the other three clusters, M62. At 
the distance of M15, the four BH candidates detected to date would 
have a radio flux density of ~ 8 Jy beam“ , which would be < 4a 
in the deepest existing image of the cluster ( [Strader et al.|2012b| ), 
and hence would remain unremarked. 

Given that both recent observational and theoretical work 
point to a population of BHs in GCs, the absence of any detected 
X-ray outbursts from such objects in the Galaxy might at first seem 
surprising. All ten of the bright, transient LMXBs in Galactic GCs 
( jBahramian et al.|2014^ have been identified as hosting NSs. This 
could be taken to suggest that either the ratio of NS to BH LMXBs 
in GCs is enhanced over that seen in the field, or the properties 
of BH LMXBs found in GCs differ from those in the field. The 
likelihood of this can be addressed empirically. Combining the last 
comprehensive LMXB catalogue [Liu et al.|2007| l with a more up- 
to-date catalogue of field BH XRBs (B. Tetarenko et al., in prep.). 


Finally, we can make a simple comparison of the relative pop¬ 
ulations of accreting neutron star and black hole X-ray binaries in 
the five clusters with deep, published radio imaging (Ml5, Ml9, 
M22, M62, 47 Tuc; [Strader et al.1|2012a|b| [Chomiuk et al.|2013[ 
this work). At present, we know of three neutron star LMXBs in 


M15 

Auriere, Le Fevre & Terzan|1984HCharles, Jones & Naylor 

1986 

Heinke, Cohn & Lugger[[2009 1 , none in M22 [Webb et al. 

2004 

Webb & Servillat[[2013), four in M62 [Pooley et al. 

2003 

Bahramian et al.||2015|(, and five in 47 Tuc ([Heinke et al. 

2003 

2005b[l. M19 has no existing Chandra or XMM data, and the deep- 


est X-ray limit to date is from ROSAT, with Lx = 3 x 10®® erg s ® 


( [Rappaport et al.|1994[ , which is not sufficiently deep to detect a 
population of quiescent accreting neutron stars. Thus we know of 
12 accreting neutron stars in our small cluster sample, as compared 
to four black hole candidates. Although the sample of clusters is 
still small, the nature of the black hole candidates is still to be con¬ 
firmed, and both populations are the subject of poorly-understood 
selection effects (not least the varying sensitivities achieved in dif¬ 
ferent clusters in the two bands), this implies a ratio in the de¬ 
tectable number of accreting black holes to neutron stars on the or¬ 
der of 1:3. Given the Galactic globular cluster population of ~ 200 
LMXBs ( [Heinke et al.|2005^ , this would then imply of order 60 
accreting black holes in clusters, which could in principle be de¬ 
tected by combining deep radio and X-ray observations, almost 
doubling the known population of black hole candidates from the 
field. 
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7 CONCLUSIONS 

We have detected persistent radio emission from the known X-ray 
source X9 in the globular cluster 47 Tuc, which had previously 
been identified as an intermediate polar. The measured radio flux 
density of 42 ± 4 /rjy beam“^ at 5.5 GHz implies a radio luminos¬ 
ity of 5.8 X 10^^ erg s”'^, which is significantly brighter than the 
persistent emission seen from any cataclysmic variable. The radio 
emission appears to be roughly constant from 2010 to 2013, as well 
as within an individual 10-hour observation. Together with the ab¬ 
sence of significant circular polarization, this also allows us to rule 
out flaring emission from a cataclysmic variable. Although transi¬ 
tional MSP systems are known to produce flat-spectrum radio emis¬ 
sion in their accreting states, they show lower radio luminosities, 
simpler X-ray spectra, and rather different X-ray variability than 
does X9. While we cannot definitively rule out this explanation, 
we find it unlikely. Combining our radio detection with archival 
Chandra X-ray data, we find that 47 Tuc X9 lies very close to 
the radio/X-ray correlation for quiescent stellar-mass black holes. 
The similarity of its optical-through-ultraviolet spectrum to those 
of dynamically-confirmed quiescent black holes further supports 
X9 being a stellar-mass black hole in quiescence, which would 
make it the fourth black hole candidate detected in a Galactic glob¬ 
ular cluster to date. While the optical magnitude implies that the 
orbital period should be < 7h (consistent with a lower-main se¬ 
quence donor star), the high X-ray luminosity of 10®® ergs“® and 
the absence of strong Ha emission suggest that the system could 
be ultracompact, with an orbital period of order 25 minutes. Differ¬ 
entiating between these two possibilities would require broadband 
optical spectroscopy. 
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